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A linear quadratic (LQ) based controller and observer concept for a semi-active full-car model is 
implemented and compared to a skyhook controller using a real vehicle in this paper. Especially an 
observer including a new modified road model is presented, which reliably filters low frequency 
disturbances induced by ascending or descending steep hills. Furthermore, methods for the parame- 
trization of the quadratic cost function are presented and it is shown in experiments for different vehicle 
masses that ride comfort and road holding can be significantly improved by using semi-active 
suspension control. In order to take nonlinear component characteristics or suspension friction into 
account a parametrized nonlinear full-car model is presented, which is used for the determination of the 
controller cost function weights. The performance evaluation is done using a 4-poster test rig as well as 
measurements conducted on a real road. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The main tasks of a vehicle suspension system are to provide a 
high level of ride comfort by isolating the chassis mass from road 
disturbances and to improve road holding by preventing the wheel 
from loosing road contact. Several studies have shown that this 
conflict can be eased by using controlled suspension systems 
instead of passive ones. Especially semi-active dampers are wide¬ 
spread in industrial application because of their performance 
increase compared to passive suspensions and their lower energy 
consumption compared to active systems (Guglielmino, Sireteanu, 
Stammers, Ghita, & Giudea, 2008). 

Methods for controlling semi-active suspension systems reach 
back to the works of Karnopp (Karnopp, 1983; Karnopp, Crosby, & 
Harwood, 1974; Karnopp & Heess, 1991) who has proposed the 
skyhook algorithm in order to improve the ride comfort. Several 
extensions of skyhook control have been suggested in the past 
decades (Ahmadian, Song, & Southward, 2004; Guo, Xu, Xing, & 
Pan, 2006; Hong, Sohn, & Hedrick, 2002; Rao & Narayanan, 2008; 
Savaresi & Spelta, 2009; Shamsi & Choupani, 2008) and especially 
the combination of a skyhook-based controller with a groundhook 
algorithm promises additional improvements on ride safety 
(Hudha, Jamaluddin, Samin, & Rahman, 2005; Valasek & Kortiim, 
2000; Yi & Song, 1999). 
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As the skyhook algorithm seems well studied, it is not surprising 
that this concept is mainly used if semi-active suspension controllers 
are implemented in real vehicles. Examples where skyhook control 
has been investigated experimentally can be found in Venhovens 
(1994) for a Volvo 480, in Lauwerys, Swevers, and Sas (2005) for an 
Audi A6, in Nieinz (2006) for an Opel Astra or in Sankaranarayanan 
et al. (2008), where a Ford Transit Connect has been tested with four 
independent skyhook/groundhook controller combinations. Also the 
control concepts of series production vehicles are based on skyhook 
control, as it can be found for Audi (Schopfel, Stingl, Schwarz, Dick, & 
Biesalski, 2007; Schwarz, Biesalski, Schopfel, & Stingl, 2008), BMW 
(Jautze, Bogner, Eggendinger, Rekewitz, & Stumm, 2008) or Mercedes- 
Benz (Friih et al., 2009). 

Even though it has been shown that clipped optimal control offers 
performance benefits for the control of semi-active suspensions 
(Butsuen, 1989; Gordon & Best, 1994; Unger, Koch, & Lohmann, 
2011 ) only quite few results on the experimental validation of model- 
based control techniques in a real vehicle are known in the literature 
so far. For the suspension control using a semi-active full-car model 
only some simulation studies exist (Jahromi & Zabihollah, 2010; Lu & 
DePoyster, 2002). 

State observers for vehicle suspension systems have also been 
investigated. In Koch, Kloiber, Pellegrini, and Lohmann (2010) a 
Kalman filter-based observer structure has been proposed which 
takes the nonlinearities of the passive damper into account. An 
experimental investigation of an fL-observer has been made in 
Aubouet, Dugard, and Sename (2010) and some studies on the 
robustness of state estimation can be found in Frohlich (2007). 

This paper extends the results presented in Unger, Schimmack, 
Lohmann, and Schwarz (2011) where the experimental evaluation of 
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an LQ-based semi-active estimator and controller concept in a real 
vehicle has been investigated. The parametrization of the quadratic 
cost function is addressed and two different controllers are presented. 
One controller is optimized to improve ride comfort considering a 
comfort measure proposed by Hennecke (1995) and the other 
controller is optimized for reducing the braking distance. Additionally 
to Unger, Schimmack, et al. (2011) a new modified road model is 
presented, which is used for the observer design in order to 
compensate for low frequency disturbances induced by ascending 
or descending steep hills. Measurements taken on a real street and 
the validation of the observer on a 4-poster test rig prove the high 
estimation quality of the concept also in challenging situations. 
Furthermore, it is shown by measurements that both controllers 
outperform the passive suspension and a skyhook-based controller in 
their respective aims, even for an increased vehicle mass. A full-car 
model, which is parametrized by measurements on a 4-poster test rig, 
is used for the observer and controller design. Additionally, an 
identified nonlinear full-car model is presented, which is used for 
the optimization routines to parametrize the controller. 

The remainder of the paper is organized as follows: In Section 2 
the nonlinear and the linearized full-car models are presented. The 
structure of the controller and the parametrization are described in 
Section 3, whereas the observer including the modified road model 
is content of Section 4. The experimental results are discussed in 
Section 5 before the paper is finally concluded in Section 6. 


2. The modeling of the system 

In order to apply model-based control techniques to a real 
vehicle, a model is required for the simulation, the controller 
and observer design and the optimization of the cost function 
weights. A nonlinear full-car model, which is presented for 
purposes of simulation and optimization, is linearized and reduced 
to seven degrees of freedom for the observer and controller 
design. The road disturbance can be described by four first order 
shape filters for each wheel which are joined with the linear full- 
car model. 

2.1. The nonlinear full-car model 

The vertical dynamics of a vehicle can be described by a full-car 
model (Chalasani, 1986), which is shown in Fig. 1. It consists of 
four wheel masses m^, each with one vertical degree of freedom, 
the chassis mass (mass m s , moments of inertia / 5 ,/„ I sS ), with 
freedoms of motion in the vertical, roll and pitch directions and 


-0.0 A 



— 1.2 A 

— 1.8 A 









relative velocity 


end stop 





deflection 


Fig. 2. The characteristics of the semi-active damper u(x,A) and the spring/ spring (x) 
(scaling omitted). 


finally the mass of the aggregate 1 (mass m agg , moments of inertia 
l a gg,/,, laggo) also with three degrees of freedom in the vertical, roll 
and pitch directions. The vehicle to be considered in this work is 
equipped with steel springs and four semi-active dampers. The 
springs (spring constant k s ) have a typical progressive character¬ 
istic including an end stop which has been directly measured and 
can be seen in Fig. 2. The dampers are of hydraulic double-tube 
type with an internal controllable solenoid valve. By applying a 
current X between 0 A and 1.8 A to the valve, the damping force 
can be modified between a minimum and a maximum bound as 
shown in Fig. 2. The tire (spring constant k us and damping rate bus ) 
can be considered as an air spring (Mitschke & Wallentowitz, 
2004) and therefore has a progressive characteristic. Moreover, 
friction effects (friction factor pi) in the suspension are considered 
using an arctangent function. All other elements are calculated 
linearly with the respective elements shown in Fig. 1. The 
equations of motion can be summarized using the suspension 
force F Sii =/ springi (x)+ /friction,i(x) (ie[1..4]), the tire force F mi , the 
forces F abr , F ab f generated by the rear and front anti-roll bars 
(stiffnesses k abr , k ab p) and the forces and torques of the aggregate 
Faggz, Tagg/p , Taggo- Furthermore, the forces of the semi-active 


1 The components engine and transmission are denoted as aggregate. 
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damper u(x, A), which are used as control inputs of the system and 
the state vector x = [z, 2:3 z 2 z 2 ■■■ z agg0 z agg0 ] T are defined. 

This leads to the dynamic equations for the wheels 


m usl z, = F usl -F sl —F ab f—li-j, ( 1 ) 

m us2 Z2 = Fus 2 -F s2 + F ab f-U 2 , ( 2 ) 

m us3 Z3 = F us3 -F s3 + F abr -U3, ( 3 ) 

mus4Z4 = F us4 -F s4 -F abr -u 4 , ( 4 ) 

the body 

m s z z = F s i 4- F s 2 4 - F s 3 + F s 4 —F agg z + Ui 1 Uj 1 u 3 1 u 4 , ( 5 ) 

4 *z* = (-F s i-Ui + Fs 2 + u 2 - 2 F abf )t / 

+(F s3 + U 3 -F s4 -u 4 - 2 F abr )t r + —T agg 0, (6) 

4 flZ 9 = (-F sl -u t -F s2 —U 2 )lf + (F s3 + U 3 + F s4 + u 4 —)! r 

+Fagg Z / ag g-T agg0 , ( 7 ) 

and the aggregate 

m agg z a gg Z = F a ggz, (8) 

i agg ,/,z a gg,/, = T agg ^, ( 9 ) 

Wagg» = W (10) 

2.2. Linearization and order reduction 


As a linear model is used for the observer and the controller 
design, the nonlinear model of the previous section is linearized 
by substituting the nonlinear elements of the model by linear 
ones. According to Mitschke and Wallentowitz (2004) and the 
parameters of the aggregate, it can be shown that the aggregate 
mainly influences the body movement at frequencies greater than 
9 Hz. In ISO (1997) it is remarked that the relevant frequency range 
for ride comfort is up to 8 Hz and therefore the degrees of freedom 
for the aggregate are neglected for the linear model. The equations 
of motion can be summarized using the disturbance vector 
a> = [Wf ■w 4 ] r and the matrices AeR (I 4xI4) , B u eR ( 14x4) and 
B 0 ,eR ( 14x4) to the state space notation 

x = Ax + B u u + B w (o. (11) 

The measured outputs y are, according to the sensor config¬ 
uration of the vehicle, the four relative suspension velocities, the 
heave body acceleration and the roll and pitch rates. The output 
equation can be calculated by using the matrices C y eR (7xl4) , 
D yu ER( 7x4 > and D yw ER< 7x4 > to 

y = C y x + Dyull + Dy w a + v (12) 

where v denotes the measurement noise, which is needed for the 
observer design described in Section 3. The parameter values of 
the linear elements are determined by optimization, as it is 
described in Section 2.4. 

For the LQ-based controller design the penalized variables are 
the vertical acceleration z z , the roll acceleration z,j, and the pitch 
acceleration z e of the body as well as the tire deflections 
Az ws>i = z—fflj by means of the scalar weightning factors p z , p„, p r/l 
and p ws . By using the matrices C w eR ( 7x14) , Dw„eR ( 7x4) and 
D w ,„eR ( 7x4) the equation for calculating the above mentioned 
outputs 


2.3. Road model 

In ISO (1995) it is reported that the power spectral density 
(PSD) of the road displacement can be approximated by 

Sw(fl=;*A00 n (14) 

with the vehicle velocity v and the roughness factor of the road A 
in [m]. As described in ISO (1995) typical roads can be classified 
using the scalar factor A, whereby higher values of A indicate 
rougher roads. The measurements in ISO (1995) show that the 
exponent n is independent from the road class, whereas the 
approximation of rt = -2 holds in general. By incorporating that 
the amplitudes of very low frequencies /-> 0 stay finite and 
considering n = - 2, Eq. (14) results in a first order shape filter 

®(t) = -/?v®(t) + v f (15) 

where (Sv is the cut-off frequency and f is a gaussian distributed 
white noise input with the spectral density A (Caspar & Nadai, 
2007; Hrovat, 1997; Koch, 2011; Venhovens, 1993). The parameters 
are choosen according to Venhovens (1993) to 0 = 0.2, v=30 m/s 
and A results from the optimization routine of the Kalman filter. 

By increasing the filter order of the shape filter, the approx¬ 
imation of the PSD can be improved especially at high frequencies, 
like it has been investigated in Turkay and Akcay (2005). But as the 
improvements are only slight and an increased filter order leads to 
higher computational effort in real-time control, the first order 
approximation of Eq. (15) is considered for the observer design in 
Section 4. 

In order to augment the system equations with the road model, 
Eq. (15) is extended to all four wheels 

a> = [~0v, -0v, -pv, -pv]U (o + [v, v, v, v]Lig (16) 

As B, 


which is then appended to the state space representation leading to 



W = CwX + DwuU + Dv, 

can be introduced. 


(13) 


Fig. 3. Experimental 
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2.4. Parameter identification 


be clipped according to (20), can be calculated as 


It is assumed that the vehicle is symmetric along its long¬ 
itudinal axis that the geometry is known and that the spring and 
damper characteristics are determined by separate component 
tests. The remaining parameters of the linear and nonlinear full- 
car models have been identified by measurements of the experi¬ 
mental vehicle on a 4-poster test rig (Fig. 3). Sinus sweeps 
between 0.5 Hz and 25 Hz for exciting the heave, pitch, roll and 
warp movements have been used as input signals for the identi¬ 
fication. As still the amount of unknown parameters is too large for 
global optimization, the parameters are combined to groups, 
which are identified separately using a Trust-Region-Reflective 
algorithm provided by the MATLAB optimization toolbox 
(Coleman & Li, 1996). The procedure is repeated iteratively until 
the variation of a global fitness value has stagnated. A comparison 
between the measured vertical acceleration and the simulation 
result is given in Section 5.2. 


3. Semi-active suspension control 

For the suspension control of the full-car, four semi-active damping 
devices are used. Accordingly, the damping force can be adjusted by a 
controllable valve, but it cannot exceed the minimum or maximum 
bounds which are shown in Fig. 2. As most semi-active control laws do 
not explicitly account for these constraints, the demanded forces u des 
must be clipped to match the controllability range of the damper, i.e. 

u = sat(u des ). (20) 


Usky=- 


WUf) “4 If ~2(! r +ff) 

2Ur+0 4E f 

2(1 !+lf) W r 2 (Ir+lf) 

’SSr^f) _4 S . 


0 0 1 ri z i 

&sky,^ 0 Z e 

0 bsky,sj 


( 21 ) 


using the geometric transformation matrix T mc - 


3.2. Clipped optimal control 

Clipped optimal control (coc) is a method where state feedback 
control can be used for the control of semi-active suspensions. 
First, a linear quadratic regulator is designed with regard to the 
fully active suspension system (a suspension system with a high- 
bandwidth actuator incorporated in parallel to the main spring 
and damper, see, e.g. Mitschke & Wallentowitz, 2004). Subse¬ 
quently, the desired forces are clipped to the controllability range 
of the semi-active damper. A state feedback law 
Ucoc = -I<cocX (22) 

is designed for the system described by (11 )-(13) which minimizes 
the cost function 

J= J o (w t (t)Q_ w w(z) + u r (Y)R w u(r)) dt (23) 

with 


Since the system is dissipative for any force generated by the semi¬ 
active damper, the overall system remains stable despite the control 
input is saturated. However, the performance may not be optimal 
for the semi-active case any more. This issue has been addressed by 
Butsuen (1989), where the optimal control law for a semi-active 
suspension system has been derived considering a quadratic 
performance index. Using the calculus of variations, Euler's equation 
and Lagrange multipliers in order to consider the state dependent 
constraints, the optimal control law for semi-active quarter-car 
suspensions has been proven to be the same as clipped optimal 
control, which is introduced in Section 3.2. For a semi-active full-car 
model a difference between the optimal control law and clipped 
optimal control can be observed, but the differences are negligible 
(Butsuen, 1989). 

Also the generalized benchmark of semi-active suspension con¬ 
trollers in Unger, Koch, et al. (2011) has been shown in experiments 
that clipped optimal control is a reasonable choice regarding compu¬ 
tational effort and performance for real-time application. Therefore, 
clipped optimal control is applied for the semi-active suspension 
control in this work. A simulative benchmark of semi-active control 
techniques considering for instance also LPV or H, control can be 
found in Poussot-Vassal, Spelta, Sename, Savaresi, and Dugard (2011). 

In order to benchmark the clipped optimal control scheme, a 
skyhook controller is used as a reference. In the following section, 
the skyhook algorithm for a full-car model is introduced, then 
clipped optimal control will be presented. 


3.1. Skyhook control 

The main idea of skyhook control, introduced by Karnopp et al. 
(1974), is to damp the sprung mass with respect to an inertial 
reference frame. In the case of a full-car model also the pitch and 
roll movements are damped inertially by applying a damping 
factor on the pitch and roll rates. The desired forces, which have to 


Qw= o 


Pa 

0 

0 


0 

0 

Pe 

0 


0 

Q , R w =p u l4. 
P K 


(24) 


The feedback matrix K coc is calculated by solving the Riccati 
equation (Anderson & Moore, 2007) 

A t P c + P C A-(P C B„ + N c )R-'(BjP c + Njj + Q. f = 0 (25) 

with P c = Pj > 0 and 

Q.c = C^CLwCw, 

Rc = DJ vu Q. w D wu + R w 

N c = c£q w D wu . (26) 

The state feedback matrix Kc 0C eR ( 4><14 > is obtained by solving 
I<coc = Rf' (B^Pc + Nj). (27) 

The force which is applied by the semi-active damper must be 
clipped analog to Eq. (20). By using the damper characteristic 
(Fig. 2) the resulting control force u is then transformed into a 
current, which is applied to the valve. 

The result is a linear control law, which is implemented in the 
experimental vehicle described by the nonlinear full-car model 
presented in Section 2.1 . The nonlinearities of the system are the 
damping force, the spring force, the tire force and friction effects. 
The main nonlinearity is the damping force, which is used as 
control input for the control design, whereby the effect of satura¬ 
tion has already been discussed. By looking at Fig. 2, it can be 
observed that the spring force, as the second important nonli¬ 
nearity, can be linearized with marginal loss of precision when the 
suspension operates within its limits. As the remaining nonlinea¬ 
rities only have minor influence on the system behavior, no further 
adjustments are made for applying the linear control law to the 
real system. 
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3.3. Determining the cost function weights 

It has been shown in Unger, Schimmack, et al. (2011) and Do, 
Sename, Dugard, and Soualmi (2011) that the parameters of the 
quadratic cost function (23) can be determined by minimizing 
more complex performance criteria using genetic algorithms. 
Thereby, two performance measures are introduced, one for ride 
comfort and one for ride safety, which are evaluated using the 
nonlinear full-car model. Then two controller parametrizations are 
optimized in order to minimize the respective measures. 

For the design of the comfort orientated controller the comfort 
index proposed by Hennecke (1995) 

ffcomf = \J fci dt) ' +k 2 Qj\l/t)dtJ ' (28) 

is used including the frequency weighted heave (a WjZ ) and pitch 
(a W/(1 ) accelerations of the body according to ISO (1997) and 
weighting factors k\, k 2 for stochastic vibrations in the time 
interval [0; t e ]. It is noted that despite the frequency range between 
4 and 8 Hz is emphasized by the weighting filter, the other 
frequencies, e.g. the range between 1 and 2 Hz which makes 
passengers feel sick, are also taken into account. 

In order to design a controller which improves the ride safety, 
the braking distance from 85 kph to a full stop on a rough road is 
used as a safety index. The results presented in Unger, Schimmack, 
et al. (2011) have shown in experiments that the braking distance 
can be reduced by 0.3 m compared to the passive hard setting and 
1.8 m compared to a lightly damped suspension setting for a 
luxury car. The procedure incorporates the extension of the non¬ 
linear full-car model by a nonlinear, longitudinal tire model 
proposed by Pacejka (2005) which describes the coherence of F z , 
the longitudinal tire force F x and the wheel slip k, i.e. 

F x = (D sin (C arctan(BX-E(BX-arctan(BX))))) + S v (29) 

with X = k + Sh, Si, = Sh(F z ), S v = S V (F Z ), D = D(F Z ) and B = B(F Z ). For 
the deceleration of the car the maximum possible force F xmaxi is 
taken, which is equal to the sum D + S v in (29) and the traveled 
distance /Z safe in [m] until the vehicle comes to a stop is evaluated. 

The tuning of the controllers is accomplished using numerical 
optimization of the unique diagonal elements of Q. w i eM (7x7) and 
R Wj ,eR ( 4x4) in (23). Thereby, in comparison with a direct numerical 
optimization of the matrix elements of Kc 0c , the number of 
optimization variables is reduced from 56 to 5 by using the 
proposed LQ-based design approach. The objective functions are 
f7comf(Q.w,i>Rw,i) for the comfort orientated controller and 
/7 safe(Q.w,2- R w,2) for the safety orientated controller respectively. 
For each function call, a closed-loop simulation of the nonlinear 
full-car model is carried out, whereby a measured country road 
profile is employed as disturbance. The optimization problem is 
solved using the genetic algorithm NSGA-II described in Deb, 
Pratap, Agarwal, and Meyarivan (2002) with the number of 
generations set to 150 and a population size of 250. The evaluation 
of the 37 500 function calls takes approximately 15 h computa¬ 
tional time on a Core i7 2.8 GHz CPU for each controller setting. 
This clearly indicates that the reduction of the optimization 
variables has been a necessary step. 


4. State observer design 

The application of the full state feedback law presented in the 
previous section to a real vehicle requires the knowledge of all 14 
states of the full-car model. As the direct measurement of all states 
is either technically hardly practicable or too expensive the 
required signals have to be estimated. The sensor configuration 
of the vehicle allows a filter-based estimation of the remaining 



states, but as it has already been investigated in Koch, Kloiber, and 
Lohmann (2010), an observer-based concept provides superior 
performance, when the estimation of the whole state vector is 
considered. Furthermore, a model-based observer allows the 
compensation of noise and disturbances on the sensors without 
introducing phase distortions. Because of the stochastical nature of 
the sensor noise and the road disturbance, the application of a 
Kalman filter algorithm is a reasonable choice. 


4.1. Kalman filter 

In order to consider the nonlinear damping force, which 
depends on the relative velocity between the chassis and the 
wheel mass and on the current applied to the valve, an observer 
structure is employed, which has been proposed by Ohsaku (2000) 
and Koch, Kloiber, Pellegrini, et al. (2010). As shown in Fig. 4, the 
nonlinear damping force is calculated externally and is considered 
as a known force input. Therefore, a linear observer structure can 
be used and the linear passive damping b si is set to 0. 

As it has been pointed out in Section 2.3 the disturbance input 
of the road model f is zero-mean and gaussian. The measurement 
noise v is acting on the output equation and is also assumed to be 
zero-mean. Both inputs are uncorrelated satisfying 

E(?? r ) = Q.{> E(^'j = R„, E(f»/') = 0. (30) 

According to Kwakernaak and Sivan (1972) the observer equation 
for the extended state space equations (17)—(19) which minimizes 
the error between the estimated state x exc and the real state x ext is 

iext = (Ae Xt -LC y , ext )x ext + (B u , e „-LD yu , c «)u + Ly (31) 

with the optimal gain 

L = P„C 7 e;(t R; 1 (32) 

and P 0 being the solution of the Riccati equation 
Ae«P„ + P„A 7 xt —PoC^R; 1 C y , ext Po + do = 0 (33) 

with Q 0 = B W ;<rQ.£Bje Xt and P„ = P„ > 0. 

4.2. Low frequency disturbance compensation 


In order to tune the Kalman filter, the diagonal elements of the 
matrices 


0. = q ( U; 


Til* 0 

r 2 0 0 

0 0 r 3 0 

0 0 r 4 


(34) 


have to be adjusted, containing the scalar weighting factors q ( , r,, r 2 , 
r 3 and r 4 . This can be either done iteratively (e.g. Frohlich, 2005; 
Venhovens, 1993) or by means of optimization (e.g. Koch, Kloiber, 
Pellegrini, et al., 2010; Unger, Schimmack, et al., 2011) using reference 
data generated by a detailed simulation model and a zero-mean road 
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Fig. 5. Absolute body velocity while ascending a road of 18% slope. 


profile as input. As already observed by Venhovens (1993) a well 
tuned Kalman filter in simulation is prone to estimate low frequency 
disturbances which occur on real streets by leaning road surfaces, 
ascending or descending hills or steady-state cornering. Fig. 5 shows 
the measured absolute body velocity while ascending a hill of 18% 
slope with the experimental vehicle shown in Fig. 3 for an observer 
laid out by the above mentioned procedure. The reference velocity 
has been generated by integration of the vertical body acceleration 
and highpass filtering with zero-phase using MATLAB's function 
filtfi.lt for eliminating the low frequency disturbance. It can be seen 
that the estimation of the Kalman filter is considerably corrupted and 
cannot be used for suspension control when no adaption is made. 

In order to get rid of the low frequency disturbance in the 
estimation, it might be possible to apply a highpass filter to the 
sensor signals or the estimated variables of the Kalman filter. The 
cut-off frequency of the filter has been selected by evaluating 
measurements taken on several different roads. The results show 
that a reliable filtering of the disturbances requires a minimal cut¬ 
off frequency of 0.8 Hz, what is in accordance to the value of 
approximatively 0.5 Hz reportet in Venhovens (1993). However, a 
cut-off frequency of 0.8 Hz causes a considerable phase distortion 
at the body's natural frequency that leads to counteracting control 
actions in semi-active suspension control. 

A large improvement can already be achieved by adapting the 
optimization routine and considering reference data measured while 
driving on a real street including hills instead of using reference 
signals generated by simulation. The reference data can be obtained 
from the internal sensors by integration and geometric considerations 
using highpass filtering with zero-phase. The optimization follows the 
procedure described in Section 3.2 for tuning the controllers. The five 
optimization variables are optimized by using the NSGA-II-algorithm 
with 150 generations and a population size of 250. The improvement 
of the estimation using the adapted optimization is shown in Fig. 5 
where it can be seen that large parts of the low frequency disturbance 
are eliminated. 

Additional benefit can be gained by assuming that the mea¬ 
surement noise is not white noise but can be modeled by 


v=H (35) 

where fi is white noise (Venhovens, 1993). This leads to a Kalman 
filter with an internal proportional-integral feedback and Venho¬ 
vens has shown that the low frequency drift in the absolute 
velocity can be eliminated for a double-lane change maneuver in 
simulation. As depicted in Fig. 5 the estimation quality is improved 



frequency in Hz 

Fig. 6. Transfer functions of the road models. 


by using the approach of Venhovens, but some disturbance effects 
still remain. 

The solution proposed in this work uses an extension of the 
road model introduced in Section 2.3 by a highpass filter with a 
cut-off frequency of f hp = 0.8 Hz. Therefore, only high frequency 
road inputs are assumed by the road model and zero mean 
estimates are enforced. In order to obtain the modified road 
model, the standard road model (15) in frequency domain 

r, - %f 

lp : fiVftg 


is extended by a highpass filter 




T hp S + 1 


(36) 


where Th P = 1 /2nf hp and ‘s’ denotes the Laplace variable. This leads 
to the modified road model 


Gstreet = C to G iD = 


(T hp S + 1)(S + pv Fz g) 


(37) 


which is applied to the full-car equations (11H13) instead of the 
standard road model. A comparison between the standard and the 
modified road model can be seen in Fig. 6 where it becomes 
apparent that both models differ in the frequency range below 
approximatively 1 Hz and are identical for higher frequencies. The 
Kalman filter is then tuned using the adapted optimization routine 
and the modified street model. 

Looking at the estimated absolute velocity of the body in Fig. 5 
it can be seen that the low frequency disturbance is eliminated 
very well and only a small error between the reference and the 
estimated velocity remains. 
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5. Experimental results 

The observer and the two controllers have been implemented 
in an Audi Q5, which is by default equipped with four relative 
displacement sensors, a body heave accelerometer, two angular 
velocity sensors for the measurement of the pitch and roll velocity 
and a stock electronic control unit (ECU). This ECU has an Infineon 
TriCore processor and is programmed using the production 
code generator dSPACE TargetLink. It is noted that no adjustments 
have been made compared to the series production semi-active 
suspension system concerning the hardware or the software 
generation process. The skyhook-based series controller has been 
replaced by the model based approach, presented in this work, 
without increasing the computational power or the amount of 
sensors. 


5.J. Observer validation 

For the observer validation the experimental vehicle has been 
excited by a road profile on the 4-poster test rig. The reference signals 
are generated using four accelerometers mounted on the comers 
of the body and one accelerometer in the middle of each wheel 
by integration and highpass filtering methods. Fig. 7 shows a 
comparison between the reference signals and the estimated states 


by the Kalman filter, 2 whereas the error between the signals is always 
very low indicating a high quality of estimation. The overall fit value 

r 2 = 1 "l4 < ? 1 ' (38) 

for this road profile is 0.65, which is similar to the values given in 
Frohlich (2007) for a real vehicle. This clearly shows that the observer 
including the modified road model operates very well in isolated 
conditions and even provides good estimates while ascending steep 
roads like it has been shown in Section 4.2. Therefore, the observer is 
used for LQ-based semi-active suspension control which is experi¬ 
mentally investigated in the next section. 

5.2. Controller validation 

The performance of the clipped optimal controller applied to 
the experimental vehicle is investigated on the 4-poster test rig 
(Fig. 3) using force sensors for measuring the dynamic tire loads 
and four accelerometers on the corners of the body for the 
measurement of the heave, pitch and roll accelerations. In order 
to evaluate the transfer behavior for the passive and the dosed- 


2 The tire positions and tire velocities front right and rear left are not displayed 
for reasons of clarity. 
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Fig. 8. Measured transfer behavior (scaling omitted). 


loop configurations, the vehicle has been excited by sinus sweeps 
between 0.5 Hz and 25 Hz. Although the closed-loop system 
is nonlinear, the investigation of the frequency response is a 
common procedure in academic and industrial research for 
the performance evaluation of semi-active suspension systems 
(Savaresi, Poussot-Vassal, Spelta, Sename, & Dugard, 2010). The 
input signals have been generated to excite the heave (all four 
hydraulic actuators move synchronously), pitch (the front actua¬ 
tors move in opposite direction to the rear ones) and roll (the left 
actuators move in opposite direction to the right ones) reaction of 
the vehicle. After a fast Fourier transform of the input and output 
signals, the transfer behavior from the road excitation at to the 
body heave (is), roll (zg) and pitch (z 7 ) accelerations as well as the 
tire force (F zdyn ) are calculated. 

The results are shown in Fig. 8 where it can be seen that both 
controllers decrease the amplitudes of the heave, pitch and roll 
accelerations over the whole frequency range. The controller 
comfort which has been optimized for ride comfort behaves like 
the softly damped passive suspension in the comfort-relevant 
range between 4 and 8 Hz (see Fig. 8, |G heaveAc cl> iCpitchAcd)■ Looking 
at the body heave and pitch natural frequencies, it can be observed 
that the peaks of the passive soft suspension are reduced by the 
controller. Therefore, the body is well stabilized by simultaneously 
providing the highest possible comfort at high frequencies by the 
comfort orientated LQ-based controller. The performance is close 
to the results presented in Savaresi et al. (2010) where the optimal 






Fig. 9. Measured and simulated body acceleration for a real road profile. 
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Fig. 10. Conflict diagrams for a road profile (left) and a sinus sweep (right). 


frequency response for a semi-active suspension system has been 
derived using model predictive control. As the body should be 
isolated from leaning road surfaces in the roll direction, the roll 
transfer behavior (see Fig. 8, I G roUAcc |) of the controller comfort 
is also satisfying because a lightly damped passive suspension 
behaves best for this usecase. 

Comparing the LQ-based controller it comfort to the also 
comfort orientated skyhook controller skyhook, the worse damp¬ 
ing properties of the skyhook algorithm in the range of the 
heave natural frequency becomes apparent (see Fig. 8, |G heave Accl. 

I GpjtchAcc I )• An increasing skyhook damping factor worsens the 
isolation properties in the comfort-relevant range, whereas lower 
values lead to an even worse behavior for lower frequencies. 
Therefore, the trade-off between body stabilization and high 
frequency comfort is tighter for the skyhook controller than for 
the LQ-based controller. 

The controller safe reduces the body movement at low frequencies 
and the tire force at the wheel hop mode considerably (see Fig. 8, 
iCheaveAccl. |G pi tchAccl. I GureForce I )• This behavior is beneficial for a 
braking maneuver because both effects have a bad influence on the 
dynamic tire load while braking. Especially the low frequency body 
movement causes huge tire load fluctuations leading to a significant 
extension of the braking distance what has been shown in Unger, 
Schimmack, et al. (2011). Considering the results of the study, a low 
rms value of the dynamic tire load stands for a low braking distance 
only if the damping of the body's natural frequency is high. As the low 
frequency damping of the controller safe is similar to the passive hard 
configuration, it cannot be improved any more. Therefore, it is 
possible to evaluate the rms value of the dynamic tire load for further 
investigations. In order to achieve this beneficial behavior with the 
controller safe, the amplitudes in the comfort-relevant range are 
increased slightly compared to the passive soft suspension (see Fig. 8, 

I GheaveAcc I < iGpjtchAccb I GeForce I )■ 

In addition to the analysis of the frequency response, which 
already contains rich information about the system behavior, the 
performance has been investigated using a real profile of a rough 
country road. The temporal responses and the power spectral 
densities of the body heave acceleration are shown in Fig. 9. In 
general, they show similar behavior as it has been observed for the 
excitation by sinusoidal disturbances. Looking at the temporal 
responses, it can be seen that neither the high amplitudes of the 
low frequency body movements of the passive soft setting nor the 
sharp acceleration peaks of the passive hard suspension can be 
observed for the controlled suspension configurations. The power 
spectral density plots indicate that the damping of the body's 
natural frequency is slightly worse for the controller skyhook than 


it is achieved by the controller comfort. The LQ-based controller 
dynamic shows a well damped behavior at the body's natural 
frequency by slightly increased amplitudes at higher frequencies 
compared to the other controllers. In summary, the results, which 
had been discussed for the frequency response, can also be 
confirmed for the real road profile. 

The comparison between the measurement and the simulation 
at X = 1.8 A illustrates the validity of the nonlinear full-car model. 
Slight differences can be observed for higher frequencies, where 
the influences of unmodeled system dynamics, like bearing 
stiffnesses, increase. 

Another performance evaluation is given by the conflict dia¬ 
grams in Fig. 10, where the Flennecke comfort index and the mean 
rms value of the four dynamic tire loads F zdyn are shown for each 
configuration. The evaluation has been made for a real measured 
road profile and the sinus sweep input in the heave direction. For 
the road profile several passive configurations have been investi¬ 
gated leading to the black curve in Fig. 10 left and it can be seen 
that both controllers achieve a minimum in their respective 
performance aims. Especially for the sweep input, where a broad 
frequency range up to 25 Hz is considered, large benefits in 
comfort and safety of the controlled suspensions can be observed. 

In order to investigate the behavior of the controllers in the 
case of a mass change of the body, the vehicle has been loaded 
with 75 kg at each front seat and 150 kg in the trunk resulting in 
an additional load of 300 kg. The results displayed in Fig. 10 for the 
controlled suspensions with additional load show a moderate 
increase in the tire load fluctuations, whereas the comfort index 
is decreased. As these effects can also be observed if any passive 
suspension is loaded the controlled suspension does not lead to 
unexpected behavior when the mass is varied during operation. 


6. Conclusion 

In this paper the experimental performance evaluation of an 
LQ-based semi-active observer and controller concept has been 
investigated. It has been shown that model-based control techniques 
can successfully be implemented in a real vehicle using series 
production hardware and offer a great performance potential. 

Especially the parametrization of the controller has been 
addressed by minimizing suitable performance measures for an 
identified nonlinear full-car model and two controllers have been 
presented to enhance ride comfort and road holding. 

Furthermore, the problem of low frequency disturbances which 
occur on real streets by leaning road surfaces, ascending or 
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descending hills or steady-state cornering causing corrupted esti¬ 
mates of the observer has been discussed. A new modified road 
model has been presented in order to compensate for these effects 
and it has been proven using measurements taken on a real road that 
the suggested method reliably filters the undesirable disturbances. 
The experimental validation of the observer on a 4-poster test rig has 
shown that all states are estimated properly. The investigation of a 
reduced sensors layout is proposed for future work. 

Moreover, the performance of the controllers has been experi¬ 
mentally investigated and ride comfort improvements have been 
shown using a measure which considers the frequency weighted, 
stochastical heave and pitch motions of the car. The road holding has 
been taken into account by minimizing the braking distance leading 
to a controller which especially minimizes the vertical tire load at the 
body natural frequencies and the wheel hop mode. Finally it has been 
shown that a modified body mass does not lead to unexpected 
behavior of the controlled semi-active suspension system. Thereby, 
LQ-based semi-active suspension control is an appropriate choice for 
the implementation in series production vehicles. 
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